Revised February 28, 2007 

(Main changes are in Sections 2 and 4.) 



Relativistic Quaternionic Wave Equation II 

Charles Schwartz0 

Department of Physics, University of California 
Berkeley, California 94720 



Abstract 

Further results are reported for the one-component quaternionic 
wave equation recently introduced. A Lagrangian is found for the 
momentum-space version of the free equation; and another, nonlocal 
in time, is found for the complete equation. Further study of multi- 
particle systems has us looking into the mathematics of tensor prod- 
ucts of Hilbert spaces. The principles of linearity and superposition 
are also clarified to good effect in advancing the quaternionic theory. 
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1 Introduction 

In a recent paper pQ, we introduced the relativistic quaternionic wave equa- 
tion, with if) = ^(x, t), 



—ip i = u ■ V-0 + mip j = (i— + k g^^ + m ^ 

which we call the free wave equation, and an extended version, 


—ip i = u • V?/> + eo^ — eu • Aip i + mi/j e ieW j, (1-2) 
at 

with external real potentials (p, A, W. For either equation we found a con- 
servation law, 

g = V.J, d,, 

where 

/9 = V>, j = {y,V*i#} (1.4) 

and {X, y} = XY + FX. 

A central detail of this study is the recognition that we must allow num- 
bers (quaternions) to multiply either on the left or on the right; so we use 
the notation (a\\b)ip = aipb. 

Here we report some further results from this continuing study, including 
two versions of a Lagrangian, further exploration of many-particle states, and 
clarification of the principles of linearity and superposition. 



2 A Lagrangian 

Formerly, we stated that we were unable to find a suitable Lagrangian for 
the quaternionic wave equation. Here are some new results. 
For the free wave equation, we go to momentum space: 

^(x, t) = , n \ 3 / 2 exp(r]u pp- x)0 Pir? (t), (2.1) 
p, v ( 27r ) 1 

where p = p/p, p — |p|, and rj = ±1. Each component 4>(t) (dropping the 
labels p, rj) satisfies the equation 

ft = #>< < 2 ' 2 > 
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where Q = ir/p + km = —Q*, \Q\ = \Jp 2 + m 2 , and we find that the following 
Lagrangian 

L = (</>-#])</>*, (2.3) 

when integrated over time, does provide a suitable action. That is, variation 
of each of the four real components of does lead unambiguously to the 
specified equation of motion (12.21) . 

This Lagrangian leads us to identify a new constant of the motion, in 
addition to 4>4>*, as follows: 

j t <W* = 0, (2.4) 

and this looks like the quaternionic version (imaginary) of the Energy. 

This is unconventional in that we have multiplied the equation (in </>) by 
its conjugate <fi* on the right rather than on the left. (Alternatively, we can 
say that we should re-identify which is the original wavefunction and which is 
its complex conjugate.) We also note that this action is imaginary, not real. 
Thus, there are actually 12 real equations resulting from the variational prin- 
ciple; and they are all consistent, implying the four real equations originally 
given. 

In fact, if we write this Lagrangian as L = iLi + jL 2 + kL 3 , any one of 
the real quantiies L\ or L 2 or L 3 taken by itself is an adequate Lagrangian to 
recover the complete equation of motion (12. 2p for the quaternionic amplitude 
0. 

Trying to find a coordinate space version of this Lagrangian, we are led 
to introduce the helicity projection operators, 



Ti, = -[l-rju- V/v 71 ^], r/ = ±l (2.5) 

7T V 7T V ' = S^rf 7T V , U • V 7T V = — TjV ~ V 2 7T^ (2.6) 

and the helicity-projected wavef unctions, ij) v = ^^(x, t). Then we can write 
the above Lagrangian as 

L = J d " x E i^jf + u ■ WV - m^kjip;. (2.7) 

This looks nice, but, of course, it is non-local in x; and I cannot extend this 
to include external potentials. See Section 4 for another attempt. 
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3 Some Quaternion Identities 

Given an arbitrary quaternion ip, there is the familiar identity = ipip*. 
Here is something else: for any imaginary quaternion numbers a and (3, we 
have the identity 

{a,fW} = {W#'} (3.1) 

involving anticommutators {, }. The easiest way to prove this is to multiply 
the left hand side by ip . . . ip* / \ip\ 2 . The result is the right hand side; but the 
left hand side is real, so it is unaffected by this operation. 

For example, the conserved current j which was previously written as 
{^, ■0*U'0} can also be written as , ipiip*} . 

Writing ip = ip\ + ip 2 , we have two further identities: 

Mi + Mi = mi + mi (3.2) 

{a, + VWi} = {P, M^l + MM- (3.3) 

In a similar vein we can generalize the conservation law (II. 3p to include a 
mixed density px% — ^1^2 + V'lV'i) "with a similar expression for the current 

Jl,2- 

A special case of the result (13. 3p . following the substitution tpi — > fiipi, is 

{a, Mi ~ Mi} = 2(^2^! - il)MD- (3.4) 

Using the above formulas, we can rewrite the momentum space Lagrangian 
of Section 2 as follows: 

Li = {-y-, L} = —(f)* i 4> + ^{Q, 4>* 1 4>}. (3.5) 

4 Another Lagrangian 

We have been able to find a Lagrangian for the extended wave equation, but 
it is non-local in time. 
Define 

i> = M,t), $ = V*(x,-*) (4.1) 
e = e (t) = +l(t > 0), or - l(i < 0) (4.2) 

and assume that the potential <p is even under t —*■ —t while A and W are 
odd. 
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Then we construct the following (real) action, 

A = Jd 4 x $ ^+{^,e^(u-V+e^}+e^u-A^-^{e ieeW k,e^}\, (4.3) 

which, upon variation of the four real components of if), gives the extended 
wave equation as follows: 

(jib 

e-^- i = u • Vip + eipi/j - eeu ■ Aip i + rrnp e ieeW j. (4.4) 

For t > this is exactly the original equation (jl.2P ; and if we now change 
t — > —t, we find the very same equation for ^>(x, —t). 

One may well wonder what is the correct choice of time t — and how 
to interpret the discontinuity which this Lagrangian posits at that point. 



5 Multi-Particle States 

In Reference [I], Section 12, we saw a particular way to construct two-or- 
more-particle wavefunctions for a generalization of our quaternionic wave 
equation; and that will be explored further below. 

In response to some comments received, I did make an attempt to ex- 
pand the established mathematics of "tensor products" , as applied to Hilbert 
spaces over an Abelian field (like ordinary complex numbers), so that it might 
accommodate quaternionic quantum theory. This effort is summarized in 
Appendix A; and it may be called a limited success. 

Let me start by repeating the previous approach used, with a slight dif- 
ference of notation, which I will explain a bit later on. Here, we will make 
frequent use of the notation for two-sided multiplication by quaternionic 
numbers and functions: (a\\b)ip = aipb. 

Write the free wave equation, Eq. (Il.ip . as 

V^=[(^\\k)-(u-V\\j)+mW = 0, (5.1) 

with our regular plane wave solution written as 

^(x, t) = tpp v (yi, t)(j) = ( exp(r/u • pp ■ x)\\exp((ir)p + km)t) )<p. (5.2) 
With x standing for the four coordinates x, t, the propagator is written as 

G + (x, x') = -(6(t - OH*) E tA^ -At- t% (5.3) 
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so that we have 

VG+ = 5\x-x'). (5.4) 

Two of these definitions, Eqs. (15.11) and (15. 3p . differ from what was given 
before by the extra quaternion k seen multiplying from the right. The reason 
for this becomes apparent when we note how the Lorentz transformation is 
applied to the wave equation: 

V^C^VC, £ = exp((^u-v||z)). (5.5) 

This is important for maintaining Lorentz covariance in some other steps, 
which involved products of these T> operators. 

If we can read Eq. (15 .4p as saying VG + = 1, can we also have G + V = 1? 
We investigate: 

G + V ^(x, t) = -J d A x\6{t - t') \\k) j^y 3 ( exp( V u ■ p p • (x - x')) 

\\exp(n PtV (t - t')) )[(^\\k) - (u- V'lij) +m] ^(x',t').(5.6) 

If we execute partial integration on the (primed) space and time derivatives, 
we do get exactly the answer ^(x, t), provided that we understand one de- 
tail of the mathematical notation here. The symbol || that stands to the 
left of the (t — t') argument in the Green function does NOT mean that 
these coordinates are to be placed to the right of the time derivative oper- 
ator which follows. That "right-multiplication" instruction applies only to 
numbers (quaternions) and not to coordinates. Our mathematical notation 
may need to be improved somewhat to make this rule transparent; but we 
shall deal with what we now have. 

Let us now see what sense we can make of the most general two-particle 
wavefunction, for the free particles, written in our "nested" way: 

*(xi, t x \ x 2 , t 2 ) = exp(r]u ■ p p • x x ) ipp :r) (x 2 , h) exp((ir]p + km)t x ) (5.7) 

^p,r,(x 2 , t 2 ) = exp(r]'u ■ p' p' • x 2 ) (f>p, v ,p', v ' exp((irj'p' + km)t 2 ). (5.8) 
PV 

Let us now study p = which we understand to be a function of two sets 
of space-time coordinates. A detailed calculation, which I will not write out 
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here, leads to the nice result 



dp 



V^j, j = {^,^*u^}. 



(5.9) 



This looks exactly like the one-particle conservation law ( 11. 3D ; but it has a 
broader interpretation now. If we identify p as a probability density, then 
this is a joint probability density that is conserved - looking at coordinates 
of particle number 1 - for any given wavefunction of particle number 2. Such 
a result in ordinary (complex) quantum theory would be rather obvious and 
would be interpreted in terms of the factorizeability of the theory for two 
non-interacting particles. 

For the quaternionic theory, this result is not trivial, especially in view of 
the previous difficulties encountered by other authors: the so-called failure 
of clustering. (See Reference [2]). 

This does not solve all such problems, however, since our starting wave- 
function (15.71) has the coordinates 1 and 2 ordered in a particular way. The 
result (I5.9p does not hold true for coordinates 2 in the derivative operations. 
However, if we first integrate over coordinate x 1; then we do recover a one- 
particle density which is conserved. 

6 More on Quaternionic Amplitude and Lin- 
earity 

In Reference p], Section 7 , we noted that the amplitude constant of a 
plane wave solution of the free equation stood in a particular place in the 
middle of the complete wavefunction, 



and not casually at the right or the left of the wavefunction as one might 
place it for convenience in ordinary (complex) quantum theory. One may 
ask whether this feature is generally true, even for the full interacting wave 
equation. 

Let us write the most general quaternionic wave equation as ^ = Hip. 
Then the general initial value problem can be resolved, formally at least, as 



ip = exp(r]u ■ p p • x) (j) exp((ir]p + km)t), 



(6.1) 



i){t) = exp(Ht) ip(0) 



(6.2) 
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where ip(0) may be specified arbitrarily and we still have a proper solution of 
the wave equation. Note that we have suppressed the space coordinates, so 
this initial data ip(0) can really be any function of space. Also, the operator H 
is likely to involve "right-multiplying" quaternions, so this ^(O) really stands 
in the midst of a complicated set of other things, just like the constant </> in 
(16.11) . If the operator H should contain time dependent terms, we also know 
how to rewrite the propagator as a time-ordered exponential of the integral 
Jo* H(t')dt'. 

This initial quaternionic data ^(0) generalizes the simple constant <fi 
which we called the amplitude of the plane wave; and from this identification 
we see the generalized meaning of linearity and superposition for quaternionic 
quantum theory. 

Let us further scrutinize the fundamental meaning of the principle of 
linearity (or superposition), which is so central to quantum theory. In the 
usual (complex) theory, we say that if ipi an d 4>2 are solutions of the wave 
equation (or state vectors in the Hilbert space), then 

i> = Ci-01 + c 2 ^ 2 , (6.3) 

for arbitrary complex numbers Ci,C2, shall also be a solution of the wave 
equation (or a state vector). When this definition of linearity is carried 
over into quaternionic theory, it causes much trouble (see Appendix A, for 
example), even aside from the question of whether to write those constants 
on the left or on the right. 

Let us ask what is really required in physics by examining the familiar 
example of the two-slit interference experiment, which is discussed in any 
textbook. Let ipi(xi) be the amplitude in the open slit number 1; and let 
^2(^2) be the amplitude in the open slit number 2. Then we want to see how 
each of those wavefunctions propagates to some point x on the observation 
screen: 

ipt(x) = G(x; £i)V>i(zi), if>2{x) = G{x\ x 2 )^2(x 2 ) (6.4) 
where the propagators G(x, derived from the pertinent wave equation. 

Finally, the principle of superposition tells us merely that we should add these 
two amplitudes at the observation screen: 

ip(x) = ipi{x) + ip 2 (x). (6.5) 

Whatever phase difference arises between these two waves (thus producing 
the interference pattern) comes from the propagators and it is thus already 
built into the correct wavefunctions at point x. 
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The mathematical statement of superposition (16. 5p is much simpler than 
the statement (16.31) ; and, I claim, this is all we need for an acceptable quan- 
tum theory. This makes a big difference in any investigation of the mathe- 
matics of quaternionic quantum theory. 
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Appendix A: Tensor Products 

A long-standing problem in earlier studies of quaternionic quantum theory 
has been the description of multi-particle states. First I will review how we do 
this in the usual (complex) quantum theory and then explore how the usual 
mathematics of tensor products may be extended to the realm of quaternions. 
We start with two particles, given coordinates X\ and x 2 - 

The usual tensor product formalism has us write a composite vector as 

ip(x h x 2 ) = f(xx)g(x 2 ) -> \i> >= f <g> g (A.l) 

This composite is supposed to reside in a new Hilbert space with the inner 
product rule 

< ^IV >=< f'\f >< g'\g > (A.2) 

which we say is factorizeable. We also note the complete linearity of this 
inner product under f C\ f\ + c 2 f 2 for arbitrary complex numbers c\, c 2 \ 
and similarly for g, f, g'. Furthermore, we have composite operators that act 
on such states as 

(A ® B)(f ® g) = (Af) ® (Bg), (A.3) 

which property we may call separability of operators. 

In the complex case, any numbers that occur as multipliers in the op- 
erators or in the state vectors can be factored out and written anyplace we 
wish, since they commute with all other operators and with one another. In 
the quaternionic case that simplicity no longer holds. In earlier study of this 
tensor product mathematics, Horwitz and Biedenharn [3] concluded that it 
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was impossible to extend this to quaternionic quantum mechanics. Their 
primary criterion was the linearity condition, which we see obviously fails 
when we have quaternionic functions / and g. This same conclusion is given 
in Adler's book [2]. 

However, as discussed above, we now have a more enlightened view of 
linearity and superposition in dealing with quaternionic wavef unctions. So 
we shall require only the simple sort of linearity: / — ► jfi + /2, which also 
allows for real multipliers. 

This still leaves us with plenty of other criteria to be checked out; and 
this is what we shall do here. 

The first place we see trouble in trying to extend the above formalism to 
quaternions is Eq. flA.20 . since the two single-particle inner products do not 
commute with one another and that violates the general rule in any Hilbert 
space 

< Tj/ty >*=< ^ > . (A.4) 



FIRST ATTEMPT 
Let us now try an alternative rule for inner products, 

< >= \[< fV ><g'\g> + < g'\g >< f\f >]. (A.5) 

which at least satisfies Eq. ( 1A.4I) . But there is more to check. Consider the 
particular composite vector 

\*o>=f®9+(fi)®(gi) (A.6) 

which, using the rule (1A.5I) . turns out to have norm < \l/ |\l/ >= 0. Yet, 
when we take the inner product < /' ® g'\^ >, we find that it is nonzero 
for general quaternionic functions / and g. This is an intolerable situation; 
and so we reject this rule flA.Sj) . 



SECOND ATTEMPT 
Alternatively, we consider the rule 



< i>'\i> >=< g'\ < f'\f > \g >, (A.7) 
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which involves one number (the inner product < f'\f >) sitting inside an 
inner product of two other vectors, < g'\ . . . \g >. For /' = /, g 1 = g this is 
simply < ip\ip >=< f\f >< g\g >. 

Searching for a zero norm composite vector we now come up with 

\*o>=f®g + {fi)®(ig). (A.8) 

Of course, for simply complex functions it is perfectly obvious that this is 
just [1 — 1]/ ® g. Is this a problem, however, for our quaternionic situation? 
We do not see the same problem as found above, in the first attempt, since, 
according to the rule (IA.T[) . the general inner product < f eg) g'\^ >= 0. 

This suggests a mathematical scheme for our tensor products where it is 
stated that / a <8> g = f <E> ag for any scalar a. This looks like the ordinary 
rule for associativity in multiplication; however, this does cause us a problem 
with the idea of separability. 

We have said that we want to have (A<S>B)(f<S>g) = (Af)®(Bg), which is 
to say that we can assign a quaternionic operator to a specific subspace. But 
we have now acknowledged that (fa) Cg> g — f Cg) (ag). These two properties 
conflict, as may be seen in 

(/ ® B) ((fa) ®g) = (fa) <g> (Bg) = f ® (aBg) (A.9) 
(I ® B)((fa) ®g) = (I®B)(f® (ag)) = f ® (Bag). (A.10) 

This is consistent only if a always commutes with B, which our quaternions 
do not admit. So this second attempt is also a failure. 



THIRD ATTEMPT 
As a variant of the second attempt, we consider the inner product rule 

< VIV> >= ~[< g'\ <f'\f>\g> + <f'\< g'\g > \f >], (A.11) 

which asserts that the order of the components /, g should not matter. Again, 
take a general superposition |\& >= / <g> g + f <g> g' and calculate the norm, 
which can be rearranged as follows. 

N 2 =< f\f >< g\g >, N' 2 =<f'\f><g'\g'>, M = NN' 
X = ±<g'\<f\f>\g>, Y = ±<f'\<g'\g>\f> 

< >= (N - N') 2 + (A.12) 
— [|1 + X| 2 + |1 + F| 2 + (1- |X| 2 ) + (1- \Y\ 2 ) . 
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Each of the five terms in Eq. (1A.12j) is non-negative, so for this norm to 
vanish, each of those five terms must vanish. The X terms in (1A.12|) will 
vanish if we choose /' = j ' \x and g' = vg with fi v = —1. This leads to the 
troublesome situation we found above in the second attempt. Now, however, 
we also need the Y terms in flA.12j) to vanish; and these two requirements 
can be met, for general quaternionic functions f,g, only with the choices 
/' = /, g' — —g or /' = — /, g' = g, which gives us trivially, |^ >= 0. 

This inner product rule (lA.lip is unfamiliar; it says, for example, that 
the vectors 1 <S> 1 and i <S> j are orthogonal to each other. Nevertheless, the 
mathematics appears to be consistent. 

So, it appears that we have now overcome some of the previous difficulties; 
and with this last choice for an inner product rule, (lA.lip . we can indeed 
achieve both the enlightened version of linearity as well as the attractive idea 
of separability given by flA.3j) . both within the confines of a proper Hilbert 
space. I do not yet have a complete proof that this is correct; but I am 
encouraged by the idea that what we required is merely a decent method of 
bookkeeping. We can keep the individual Hilbert spaces effectively separated 
but still have a proper composite Hilbert space with this tensor product. 

Thus, if we have two operators, A acting on the vector / and B acting 
on the vector g, then we can write them, in the tensor product notation, as 
A® I and / <8> B, respectively; and these two operators now commute with 
each other regardless of any quaternions embedded in them. 

How to extend this rule from two to n components in the tensor product 
is an open question. It might stay at two terms in the inner product: a 
given order and its reverse. Or it might become 2 n ~ 1 terms by a process of 
induction; or it might become n!, with all permutations of ordering. 

What about the criterion of factorizeabilityl We have clearly had to give 
up on that. I think this is more a matter of physics than of mathematics. If 
this quaternionic theory has any relation to physical reality, then it appears 
that we shall have to acknowledge some sort of persistent entanglement of 
many-particle states. That is not nice, according to our conventional ideas 
and experience; but it is not unimaginable. The example of two-particle 
density studied in Section 5, above, gives us a mixed picture of such entan- 
glement. 
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